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LOCAL RADIATION MHD INSTABILITIES IN MAGNETICALLY STRATIFIED MEDIA 
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ABSTRACT 

We study local radiation magnetohydrodynamic instabilities in static, optically thick, vertically 
stratified media with constant flux mean opacity. We include the effects of vertical gradients in 
a horizontal background magnetic field. Assuming rapid radiative diffusion, we use the zero gas 
pressure limit as an entry point for investigating the coupling between the photon bubble instability 
and the Parker instability. Apart from factors that depend on wavenumber orientation, the Parker 
instability exists for wavelengths longer than a characteristic wavelength Atrani while photon bubbles 
exist for wavelengths shorter than At ran - The growth rate in the Parker regime is independent of 
the orientation of the horizontal component of the wavenumber when radiative diffusion is rapid, but 
the range of Parker-like wavenumbers is extended if there exists strong horizontal shear between field 
lines (i.e. horizontal wavenumber perpendicular to the magnetic field). Finite gas pressure introduces 
an additional short wavelength limit to the Parker-like behavior, and also limits the growth rate 
of the photon bubble instability to a constant value at short wavelengths. We also consider the 
effects of differential rotation with accretion disk applications in mind. Our results may explain why 
photon bubbles have not yet been observed in recent stratified shearing box accretion disk simulations. 
Photon bubbles may physically exist in simulations with high radiation to gas pressure ratios, but 
higher spatial resolution will be needed to resolve the asymptotically growing unstable wavelengths. 
Subject headings: accretion, accretion disks — instabilities — MHD — radiative transfer 
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1. INTRODUCTION 



Considerable attention has been devoted recently to the study of radiation pressure driven instabilities in magnetized 
media, which have generally come to be known as ph oton bubble i nstabilities. This name is particularly appr opriate to 
the bu oyant long wavelength instabilities studied by lArons! (|1992l ) in the context of accreting X-ray pulsars. iGammie 



2001 



1998) foun d that an analogous instability is likely to exist in radiation-dominated accretion disks. iBlaes fc Socrates 



2003) extended Gammie's analysis to finite gas sound speeds and short wavelengths to show that the instability 



in this regime is a radiative amplification of magnetosonic modes. At le ast in med i a whe re Thomson scattering is 
the dominant opacity, the slow mode always has the fastest growth rate. IBegelmanl (|2001f ) derived a fully nonlinear 
periodic shock train solution to the equations of radiation MHD, which he suggested would be the nonli near outcome 
of th e short wave l ength p hoton bubble instability. This was later confirmed by numerical simulation (| Turner et al.l 
2005). Begelman (2006a) has also discovered a nonlinear wave solution in the long wavelength, buoyancy regime of 
the photon bubble instability. 

Application s of the short wavelength photon bubble instability to accretion disks have been explored extensively 
bv IBegelmanl (|2002L [2006b). Perhaps the most significant implication is that the instability might permit highly 
super-Eddington luminosities from accretion disks that would still be geometrically thin in terms of their vertical mass 
distribution. The physics of the short wavelength photon bubble instability has also been extended to situations where 
other forms of diffusive energy transport ex ist. In particular, work has been done on "neutrino bubble instabilities" 
in proto-neutron stars ()Socrates et al.|[2005l) and " Coulomb bubble instabilities" in systems with anisotropic thermal 
conduction due to the presence of a mag netic field ([Socrates. Parrish fc Sto"nell2008D . 

All of this work on photon bubble and related instabilities has assumed a uniform background magnetic field, but 
gradients in the magnetic field can also drive in stabilities. In parti c ular, a medium may be vulnerable to magnetic 
interchange and undulatory Parker instabilities (|Tserkovnikovlll960l : lNewcombi ri961: Parker 1966. fl967t ) if magnetic 
pressure gradients contribute significant support against gravity. Vertically stratified shearing box simulations indicate 
that s uch instabilities, n ot photon bubbles, appear to dominate the large scale dynamics of the surface layers of accretion 
disks (Bl aes et al1l2007l ). though it is conceivable that photon bubbles would be relevant on smaller length scales that 
are unresolved in the existing simulations. 

If both the equilibrium density and magnetic pressure decrease outward, then both the Parker and photon bubble 
unstable modes reduce to slow magnetosonic modes in the short wavelength limit. It is therefore likely that these 
two instabilities have nontrivial couplings in equilibria with both magnetic and radiation pressure gradients, and it is 
important to study this coupling in order to understand under what regimes and length scales each instability is likely 
to operate. This is the purpose of the present paper: we investigate local radiation MHD instabilities in equilibria 
that are supported against gravity by both radiation and magnetic pressure gradients. 

We begin in section 2 by stating the basic radiation MHD equations and assumptions that we employ in our 
analysis. The perturbation equations are very complicated to analyze when we include both radiation and magnetic 
field physics, and much of the rest of the paper is devoted to various levels of approximation. The situation is simplest 
when one neglects gas pressure completely, as in Gammie's (1998) original analysis of the photon bubble instability. 
We employ that approach in section 3, and find a single unstable mode at all wavelengths. There is a characteristic 
transition wavelength Atran above which the mode is Parker-like and below which the mode becomes the standard 
photon bubble instability, although the actual transition wavelength depends on the orientation of the wave vector of 
the perturbations. In section 4, we discuss the effects of finite gas pressure. Because much of the recent interest in 
photon bubble instabilities lies in accretion disks, we include the effects of differential rotation in section 5 and then 
discuss applications of our results to accretion disks in section 6. We finish by summarizing our findings in section 7. 
In the appendices we briefly show that our perturbation equations recover the basic properties of magnetic buoyancy 
instabilities derived by previous authors and show that, in contrast to the case of adiabatic perturbations (Ncwcomb 
1196 lh , the most rapidly growing Parker modes do not require large horizontal shear between field lines when radiation 
diffusion is rapid. We also present a more detailed derivation of our instability analysis with differential rotation. 

2. EQUATIONS 

IBlaes fe Soc rates (2003) derived simplified versions of the general radiation MHD equations o flStone. Mihalas fc Normanl 
(|1992t ). and we continue to use the same equations here, but with some additional assumptions. We assume that the 
gas and radiation temperatures are locked together and that the flux mean opacity k is a constant (as would be true 
for Thomson scattering). The gravitational field g = — <?z, with g > possibly being a function of height z, but 
independent of time as we neglect the effects of self-gravity. The equations then become 

^+V-(pv)=0, (1) 

p (| + v • Vv) = -V P + PS + 1* • VB - i-V^ + f F, (2) 

d /4 \ 

(e + E) + v • V(e + E) + [ -E + 7 e V • v = -V • F, (3) 



F =-£p- VE > (4) 



Radiation MHD Instabilities 



3 



_ = Vx(vxB), 



(5) 



and 



P 



P = 



pk B T 



and E — aT 



(6) 



7-1 ii 

Here p is the fluid mass density, v is the fluid velocity, p is the gas pressure, e is the gas internal energy density, T is 
the temperature, E is the radiation energy density, F is the radiation flux, and B is the magnetic field. Other symbols 
have their usual meanings: c is the speed of light, a is the radiation density constant, fee is Boltzmann's constant, \x 
is the mean particle mass in the gas, and 7 is the ratio of specific heats in the gas. 

2.1. Equilibrium 

We assume a static, vertically stratified, horizontally homogeneous equilibrium with purely horizontal magnetic field 
B = B(z)y. The equilibrium radiation flux is F = Fi. The only differential equations that the equilibrium must 
satisfy are those of hydrostatic equilibrium, 

E 

3 + 8tT 



d_ 

dz 



= -pg, 



radiative equilibrium, 



and radiative diffusion 



dz ' 



p= __^dE 
2>np dz 

2.2. Perturbations 



(7) 

(8) 
(9) 



We linearize equations (H}-© about the general equilibrium, assuming a (t, x, ^-dependence of exp[i(k x x+k y y— wt)]. 
We do not make any assumptions about the z-dependence of the perturbations as our background is a function of z 
and we in particular wish to include the effects of the background magnetic gradient to the maximum possible extent. 
Eliminating magnetic and radiation flux perturbations, we obtain five coupled ordinary differential equations: 



iujSp + ik x 8v x + ik y 8v y 



d5v z Sv z 



dz 
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dz 
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iv A 



d 2 5v z 
udSP 



and 



p dz 

AcE d 2 ST 
Sup dz 2 
^dSp 
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ujgSp = 0, 



+ i 



k 2 v 2 A + 



kyV A 
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(12) 
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Sv x 



mag 



(13) 



- 8F- 
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ZnpHp 
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P^ + ^T(cf-ci)Sp- 
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1 
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AE 

IT 



ST 



-Sv z = 0. 



dz ' r 3 -i v_I " x '" r s(r 3 -i) 

Here v A — B/^Airp is the Alfven speed, Sp = 5p/p, ST = ST/T, and SP is the total thermal pressure perturbation, 



(14) 



8P = Sp 



5_E_ 
"3" 



PCi Sp 



P 



AE\ 



ST. 



(15) 



The quantities H p (z), _ff mag (z), and H!^ (z) are measures of the local density and magnetic pressure scale heights, 



H„ 



dhip 

dz 



^mag — 



rflnB 2 

dz 



H = 

mag 
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2 dz 2 



1 /dlnB 2 

2 V dz 



-1/2 



(16) 
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We implicitly assume that H p and H mag are non-negative, i.e. that density and magnetic field do not increase outward. 
This is the regime of most relevance for the outer layers of stellar envelopes and the uppermost layers of accretion 
disks. The quantity ./V 2 is the square of the local Brunt- Vaisala frequency in the gas-radiation mixture, 



N 



9 



J__d_ 

Ti dz 



In 



P 



din p 



dz 



(17) 



c\ = (p/p) 1 ^ 2 is the isothermal sound speed in the gas, c t = \Pi(p + E/3)/p] 1 ^ 2 is t he total adiabatic sou nd speed in 
the gas plus radiation mixture, and Ti and T3 are generalized adiabatic exponents (Chandrasckhail ll967f) . 



Ti = 



16E 2 + 60(7 - l)£e + 97(7 - l)e 2 
9(e + 4E)(p + E/3) 



and I\ = 



16£; + 37e 
3(e + 4£) ' 



(18) 



3. THE ZERO GAS PRESSURE LIMIT 

The perturbation equations presented above are complicated and cannot be readily combined into a single ordinary 
differential equa tion. To make progress, we first consider short wavelengths where radiative diffusion is rapid, and follow 
Gammic (1998) by adopting the limit of zero gas pressure. After all, we are particularly interested in environments 
where magnetic and radiation pressures dominate greatly over gas pressure. The mathematical advantage of this 
approximation is that there is then no slow magnetosonic wave limit and the photon bubble instability appears at 
lowest order in the short wavelength limit. (For nonzero gas pressure, this analysis would only be valid for photon 
bubble wavelengths longer than a finite turnover wavelength; see section [4] below.) At the same time, the minimum 
unstable wavelength for the Parker instability also vanishes in the rapid diffusion limit if the gas pressure is zero. 
Hence the zero gas pressure limit enables us to explore the coupled photon bubble and Parker problem in the short 
wavelength WKB limit at lowest order. We shall find that the shortest wavelengths are always in the photon bubble 
regime, but if magnetic support dominates radiation support in the equilibrium, then there is a WKB transition 
wavelength longward of which the Parker instability becomes manifest. 

Proceeding with a WKB ansatz that the z-dependence of all perturbations is exp(ifc z z), we obtain a seventh order 
dispersion relation. Taking the infinitely high wavenumber limit with co ~ k 2 gives the usual damped diffusion mode, 



ck 2 

u = -i; — , 
6np 

where k 2 = k 2 + k 2 + k 2 . Doing the same thing with uj ~ k gives the two fast magnetosonic modes, 



(19) 



2 ; 2 2 

LU = k v A , 



and the two Alfven modes, 



,2 2 

KyV A . 



(20) 



(21) 



The remaining two modes describe the coupled Parker and photon bubble instabilities. Assuming that uj grows no 
faster than fc 1 / 2 as k — > 00, the lowest order terms in the dispersion relation become 



to 



2 k 4 



/4kE 
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k 2 y k 2 
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mag 
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Clearly as k — > 00, the lowest order solution is given by 



k 2 k 2 



2k z k y/ 



kF 



2H n 



ik z kyk 



nF 



0. 



(22) 



9 . k v k z kF 
w = i-^rr. , 



(23) 



identical to the photon bubble dispersion relation in this regime (IGammielll998f) . 

On the other hand, if magnetic pressure gradients dominate radiation pressure gradients in supporting the medium 
against gravity, the last term in equation (|22l) can be small except at extremely short wavelengths. Neglecting this 
term, and also the linear term in uj for the moment, the dispersion relation gives for large magnetic pressure gradients 



(k 2 x + k 2 y )v 2 A (k 2 x + k 2 y )g 



Ak 2 H 2 

11 mag 



2k 2 H„ 



(24) 



which is the short wavelength Parker instability growth rate in this reg ime. This result a grees with the c\ — and 
k — > 00 (with lu ~ k°) limit of the Parker dispersion relation derived bv iBlaes et al.1 (|2007h in their appendix. Their 
result assumed an isothermal medium, implying F = 0. However, we show in Appendix |A] that it is also valid for more 
general equilibria in which radiation and magnetic pressures dominate gas pressure, and radiative diffusion is rapid. 
[See equation (|A15[) .] 

Note that the pure Parker instability is a genuine exponentially growing instability in all space, whereas the pure 
photon bubble instability is an overstability of traveling waves. Nevertheless, the different wavenumber dependencies 
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of equations (|2"3")l and (|2~i)) indicate that a transition between photon bubble and Parker behavior will occur at a 
characteristic vertical wavenumber fctran given by 



2tt 



(25) 



provided at least that fc^ is not too large compared to fc y . Note that the photon bubble growth rate always declines 
with increasing k x at fixed k. The photon bubble instability grows fastest when the perturbation gradients lie entirely 
in the plane defined by gravity and the magnetic field. The Parker growth rate is relatively unaffected by k x in this 
zero gas pressure regime. 

The corresponding transition wavelength can also be written in a more physically transparent fashion as 



A-kH, 



mag 



V dP mag /dz 



(26) 



where P rac i = E/3 is the radiation pressure and P mag = B 2 /(8tt) is the magnetic pressure. As magnetic pressure 
gradient forces become stronger compared to radiation pressure forces, At ra n decreases, extending the long-wavelength 
range of the Parker instability at the expense of the short-wavelength range of the photon bubble instability. 

For WKB to make sense with k < fctran, we require fctran^mag 2> 1- This implies that magnetic support must 
dominate radiation pressure support in the equilibrium. If this is not true, then there is no WKB-accessible range 
of wavelengths where the Parker instability will be present (except for very short radial wavelengths: k 2 /k 2 S> 

(&tran-ffmag) )• Assuming fctran-ffmag ^ lj the middle term in the square brackets in equation (|2"2")l is negligible, 
leaving us with 



u 2 + i v 



klAEn {k 2 x + k 2 y )v\ k 2 y k z nF 



0. 



k 2 3c ±H 2 k 2 k 2 , 

mag 

If damping (the linear term in ui) is negligible, the solution to this dispersion relation is 

k 2 ^ 1/2 

fx ~, 



(27) 



kyV A 



(1 + *WH 

U2 \ 1/2 



k 2 

"-trail 



(i + K/K) + 



k 2 

"'tran 



(i + K/k 2 ) 

1/2 v 



1/2 



+ (1 + k 2 Jk 2 ) 



(28) 



The positive root is unstable: photon bubbles are recovered for k z S> (1 + k 2 /k 2 )kt ra _ n , and Parker is recovered for 



k z < (I + k 2 x /k 2 y )k tran . 

To understand the role of the damping term, it is helpful to introduce the dimensionless quantity 

AEnHr, 

— = 3 | 

c . 



3cva 



'(7)sM 



(29) 



where c r = (4E ' j'dp) 1 ! 2 is the effective radiation sound speed. It is also convenient to scale the angular frequency of 
the mode with fA/^mag, 

UJ = -. 



VA 

Then the solution of the dispersion relation (l27l) may be written as 



. y> ky_ 

2k 2 2k 



- 1 



h y 



k 2 



1/2 



(30) 



(31) 



or 



% 2k 2 2 3 / 2 fc' 



(1 + kl/k 2 + k 2 b 2 /k 2 ) 2 + 



k 2 

^tran 



1/2 
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1/2 
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(32) 



Note that the mode with the upper (plus) sign is always unstable, in that the imaginary part of the frequency is 
always positive. The other root is always damped. The unstable mode again transitions smoothly from Parker for 
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Fig. 1. — Scaled unstable mode growth rate 7 = Im(tj) as a function of the dimensionless wavenumber fc/fctran- The different colors are 
for different wavenumber orientations: k = (0,0.707,0.707) (black) and k = (0.9,0.1,0.424) (red). Solid curves have no damping (b = 0), 
while dashed curves have b = 10. 



k z <C (1 + kl/ky + kyb 2 / k 2 )k tra , n to photon bubble for k z > (l + k 2 /k 2 + k 2 b 2 / k 2 )k tra , n . The instability always becomes 
Parker-like for wave vectors in the x — z plane, i.e. k y — > 0. 



Damping affects the Parker regime most. For k z <C (1 + k x jk + k b < jk )fc tran , we have 



ui 



l 2k 2 



±i- 



kp 2 



1/2 



2k 



(33) 



Unless k y —> 0, the growth rate of the unstable root is greatly reduced when damping is large (d > 1): G) ~ 
i(l + k 2 /k 2 )/(4b). On the other hand, photon bubbles are relatively unaffected by damping provided fc z //ctran ^ b 2 . If 
instead b 2 3> k z /k tran ^> 1, then the real and imaginary parts of the unstable photon bubble frequency are markedly 
reduced: 

4Ek ( k z , k 2 z k 2 \ 

u ~ 3c \Ak^ + l mkl an k 2 y b*J ■ {3 ) 

Figure [1] illustrates these behaviors. The black curves show cases with k x — and zero damping (b = 0, solid 
curve) and finite damping (b = 10, dashed). Here the mode transitions from constant (Parker) growth rate at low 
wavenumbers to a (photon bubble) growth rate that increases as k 1 ! 2 for high wavenumbers. The transition occurs 
at k ~ fctran with zero damping, and for k somewhere between 10 and 100 fc tran for b = 10. The latter is consistent 
with the transition wavenumber being given by k z ~ k 2 b 2 kt ran /k 2 , implying k ~ (ky/k)b 2 k tran = 70.7fct ra n- For this 
wavenumber orientation, damping affects the Parker instability severely. The photon bubble growth rate is much less 
affected, however: little damping occurs for k z > b 2 k tra ,n or k > (k/k z )b 2 kt ra ,n = 140fctran- The red pair of curves 
have the same damping parameters but now have wavenumbers oriented mostly in the x — z plane. The transition to 
photon bubble behavior now occurs above a much higher wavenumber k ~ (k 2 J k 2 )(k j k z )k tTa , n ~ 200fct ra n- The Parker 
instability is insensitive to damping for this wavenumber orientation. 

4. FINITE GAS PRESSURE EFFECTS 

We now consider the effects of nonzero gas pressure on the coupled photon bubble/Parker instability. Finite gas 
pressure stabilizes the Parker instability at infinitely short wavelengths. Then radiative diffusion is so fast that ST ~ 0, 
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and instability only occurs for wavenumbers satisfying 



< 



~; kv 



(35) 



where 



'mag*'; 



1/2 



(36) 



[See equation (|A12[) in the Appendix.] If the magnetic pressure gradient is the dominant force in supporting the 
medium against gravity, then kp reduces to « g/(ciVA) = va/ (2H ma _gCi). In this case equation (|35p can be written as 
kyVA < g/ci (ignoring wavenumber factors), which reveals its physical origin. Buoyancy is fundamentally what drives 
the Parker instability. Stability therefore requires that the magnetic tension not have time to restore the perturbed fluid 
elements on a buoyant rise time in the gas in the ambient gravitational field (which is largely determined by magnetic 
pressure gradients). Note, however, that the growth rate of the fastest mode of the instability is ~ t>A/-ff m ag, and 
occurs at a wavenumber much smaller than the marginally stable wavenumber kp [equation (|A16[) in the Appendix] . 

In a medium that is only supported against gravity by gas and radiation pressure (no magnetic field gradients), 
finite sound speed affects the photon bubble instability by causing the growth rate to reach a fixed value once the 
wavenumber is high enough for the mode to become a slow magnetosonic mode. Provided that gas and radiation 
temperatures are thermally locked (as we are assuming throughout this paper), this turnover wavenumber is given by 



kF 



3p 
AE 



(37) 



For wavenumbers greater than /jt, the instability growth rate asymptotes to a constant value ~ fcxCi when the magnetic 
pressure is much larger than the gas pressure p (jBlaes k, Socratesll2003D . 

Previous studies of the short wavelength photon b ubble instability assumed a uniform equilibri um magnetic field, 
which therefore provides no support against gravity (jGammiel 119981 : IBlaes fc Soc rates 2001, 2003). In this case, the 
turnover wavenumber fcx is approximately g/c?, i.e. the reciprocal of the gas pressure scale height (jBlaes fc Socrates! 
2003). When magnetic pressure gradients support the equilibrium, g/cf can be significantly larger than the true 
turnover wavenumber, implying that the asymptotic short wavelengths of the photon bubble instability (which have 
the highest growth rates) can be more easily resolved in numerical simulations than one might have expected from 
just the gas pressure scale height. On the other hand, the asymptotic growth rate can also be significantly smaller 
than the radiation pressure supported equilibrium estimate g/ci when magnetic gradients support the equilibrium. 

We can write the three wavenumber scales in the problem as 
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(38) 



(39) 



(40) 



where 77 rac i = — {dlnE/dz)^ 1 is the radiation pressure scale height, &t and fcp are the photon bubble turnover and 
Parker cutoff wavenumbers, respectively. We emphasize that these expressions are only valid in the limit of rapid 
radiative diffusion. To study our problem around these wavenumber regimes with WKB methods, we need kH 3> 1, 
where k is any of fctran, &t or fcp and H is any of the scale heights. 

As we showed in section [31 our a = WKB dispersion relation (|2"T)l agrees with previously published photon bubble 
and Parker instability growth rates. This suggests that taking the z-dependence of all perturbations to be exp(ifc 2 2;) 
without combining the perturbation equations first may prove useful even in the finite but small (compared to magnetic 
and radiation pressures) gas pressure regime. We solved the resulting dispersion relation numerically. Figure [2] shows 
results for an illustrative set of parameters. The photon bubble growth rate asymptotes to a finite value beyond the 
turnover wavenumber fcx because the gas pressure is no longer exactly zero. 

In both Figures [T] and [U note that there is no sign of photon bubbles for low wavenumber (fc -C fctran) where diffusion 
will eventual ly be slow. In particular, our low wavenumber growth rate agrees with that of the Parker instability, but 
not with the lAronsI ([19921 ) photon bubble result. Thus if magnetic gradients are large enough that fctran corresponds 
to wavelengths in the rapid diffusion regime, then the slow diffusion version of the photon bubble instability may not 
exist. However, a non-WKB analysis is necessary to verify this conclusively. 

It is also possible in principle for the maximum Parker growth rate to be higher t han that of the photon bubble 
instability. Analytically, this occurs when the photon bubble gro wth r ate obtained by IBlaes fc Socrates! (|2003h (th en- 
equation 93) exceeds the Parker growth rate stated in equation (IA13|) . In the limit that the medium is magnetically 
supported, and the gas pressure is much smaller than radiation and magnetic pressures, this comparison means that 
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Fig. 2. — Instability growth rate 7 = Im(cj) (black curve) in units of fctran^A versus total wavenumber in units of fctran, f° r a wavenumber 
orientation k x : k y : k z = sin(7r/12) : 1 : 1. Here we set all the scale heights in the calculation to be equal: -ff mag = H p = -ff ra d- The sound 
speeds, Alfven speed, and damping parameter satisfy Cj/ua = 0.01, c T /v\ = 0.1, and b = 0.01. This implies ktr&nH ~ 33, fcp/ fctran — 50, 
and fex/fctran — 75. The blue curve in this figure shows the growth rate from the analytic zero gas sound speed expression l|32[ l, 

Parker dominates if VA/H mag ^> (c r /7? ra d)(c r /ci). Provided the scale heights are comparable, this is equivalent to 
the condition that kp < fctran, so that the Parker growth rate cuts off and decreases down to the asymptotic photon 
bubble growth rate. In this case the Parker instability transitions into the photon bubble instability when fc exceeds 
fcp, which is less than the transition wavenumber fctran for zero gas pressure. For finite gas pressure, the Parker growth 
rate can exceed the photon bubble growth rate even without this ordering of the characteristic wavenumbers. We will 
later present a finite gas pressure example of this from a numerical simulation of accretion disk vertical structure in 
section E] below. 

5. EFFECTS OF ROTATION AND SHEAR 

Motivated by the fact that both Parker and photon bubble instabilities are driven by vertical gradients, and by 
possible applications to stellar envelopes and atmospheres, our analysis so far has assumed a static equilibrium. 
However, there has been considerable recent interest in photon bubble physics in radiation dominated accretion disks. 
Such flows are differentially rotating, and before one can apply our resul t s to s uch flows, one has to account for 
the effects of rotation a nd shear. Such effects were first examined by iShul (|1974f ) for the Parker instability, and by 
iBlaes fc Socra tes (200ll) in the case of the photon bubble instability restricted to axisymmetric perturbations. 

In the spirit of the short wavelength WKB approximation that we have been using, we restrict attention to 
a local, corotating patch of a differentially rotating disk and u se the fluid equations describing a shearing box 
(jGoldreich fc Lvnden-Bell| [l965; Ha wlev. Gammie fc Balbus~lll995) with coordinate axes (x,y,z) aligned with the ra- 
dial, azimuthal, and vertical directions, respectively. Appendix [B] presents a detailed analysis of the coupled Parker- 
photon bubble problem, again assuming negligible gas pressure as we did in section |3] above. Provided kt ran H mag 1, 
the azimuthal component £ y of the small amplitude Lagrangian displacement of perturbations satisfies the differential 
equation 



d 2 i y k 2 y AEKdZ y (*x + *>A. . AkF 



dt 2 k 2 3c dt 



AH 2 k 2 



Cy = 0. 



(41) 



In the absence of differential rotation, all the coefficients of this equation would be independent of time. We could then 
assume a complex exponential time-dependence £ y cx exp(— iut) and recover the static equilibrium dispersion relation 
(|27[) that we derived above. As we discuss in more detail in Appendix [B] Coriolis forces are negligible provided the 
WKB approximation is valid. 

However, the shear due to the differential rotation of the equilibrium flow has a more important effect. This shear 
means that the radial wavenumber k x depends on time: k x (t) — k x g + qCltky, where k x g is the initial value of k x , q 



Radiation MHD Instabilities 



9 




Fig. 3. — Upper plot: modulus (solid curve) and real part (dashed curve) of the azimuthal component of the Lagrangian displacement, in 
arbitrary units, as a function of time in a shearing box with negligible gas pressure. This particular case assumed dimensionless parameter 
values of b = (i.e. negligible radiative damping), f!H ma g/iiA = 1, q = 3/2, fcz/fctran = 90, and (k x Q, k y , k z ) = (—0.993, 0.05, 0.107)£:. 
Lower plot: the dimensionless ratio fcyfc z /(fc^fctran) as a function of time for the same parameters. The instability is in the Parker regime 
when this ratio is less than of order unity, and in the photon bubble regime otherwise. 

is the shear parameter (3/2 for a Keplerian accretion disk), and VI is the angular velocity of the local patch of disk. 
The exponential growth rates of the static equilibrium analysis will therefore only be valid so long as these growth 
rates are faster than the time it takes the radial wavenumber to change significantly in time. For a magnetically 
supported medium, this implies that {k 2 /k 2 ){k 2 _/k 2 ) ^> q 2 for the Parker regime, and {k 2 /k 2 )(\k z \/k tran ) ^> q 2 for 
the photon bubble regime. The latter is easier to satisfy than the former, as the transition from Parker to photon bubble 
occurs when \k z \/kt Yaa > k^/k 2 , assuming negligible damping. On the other hand, if radiation pressure supports the 
equilibrium, then shear will not affect the exponential growth of photon bubbles provided {k 2 /k 2 )\k z \H va< \ ^> q 2 , 
which is essentially equivalent to the WKB condition. Short wavelength photon bubbles are therefore generally less 
affected by shear than Parker modes. 

Note that if these conditions are not satisfied and/or the instabilities fail to reach the nonlinear regime before the 
waves start to shear, then the character of the unstable driving can change with time. For example, a strongly leading 
nonaxisymmetric perturbation (e.g. k y > and k X Q large and negative) could start in the Parker regime, then flip 
to the photon bubble regime as the wave swings from leading to trailing (k x ~ 0), and then flip back to the Parker 
regime as the perturbation becomes strongly trailing (k x large and positive). We illustrate this linear regime behavior 
in Figure [3] which depicts a numerical solution to equation (|41[) . Nonaxisymmetric photon bubbles inevitably enter the 
Parker regime unless they become nonlinear first as the wavevector shears to be strongly trailing. But again, photon 
bubbles whose vertical wavelengths are short enough should be able to grow much faster than the shear rate. 

Axisymmetric (k y — 0) perturbations are immune to these shear effects, but the photon bubble growth rate vanishes 
for such perturbations. Photon bubbl es necessarily require w ave vectors which are neither entirely perpendicular or 
along the equilibrium magnetic field (jBlaes fc Socrate s 2003). Because we are considering equilibria with vertical 
magnetic gradients, we restricted consideration to equilibrium magnetic fields that are entirely in the azimuthal (y) 
direction. These are generally the dominant magnetic field components even in a turbulent accretion disk because of 
the differential ro tation. But other components are present, especially in the surface layers if large-scale Parker modes 
are able to grow (Bl aes et al1l2007t ). Axisymmetric photon bubbles on an equilibrium with magnetic field components 
in all directions have growth rates that are proportional to the square of the ratio of the non-azimuthal field components 
to the total magnetic field ()Gammiei 119981 : iBlaes &; Socrates! 120011 12003T ). and hence have slower growth rates than 
nonaxisymmetric photon bubbles in a medium dominated by azimuthal fields. 

Our discussion in this section has assumed negligible sound speed c; in the gas alone. As we will see in the next section, 
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radiation MHD simulations of local patches of accretion disks generally produce conditions in which fcx > fcp > fctran- 
Hence the main thing that finite gas sound speed introduces is a limit ~ fcxCi to the short wavelength photon bubble 
growth rate for k J> /ct- In that very short wavelength regime, photon bubbles continue to be immune to the shear 
rate. For a radiation pressure supported medium, this merely requires Q,H Taj ±/c\ 3> q, which is clearly satisfied 
because ~ c r /7J ra d. For a magnetically supported medium, this requires instead that ^/(qua) 3> 5 provided the 
radiation pressure and magnetic pressure scale heights are comparable. This in turn is equivalent to the condition that 
/ct 3> fctran, which we find to be generally true. 

To summarize, in a medium supported by radiation pressure, rotation and shear do not significantly affect the photon 
bubble regime provided we are considering vertical wavelengths that are short enough for the WKB approximation to 
be valid (\k z \H ^> 1). In a magnetically supported medium, photon bubbles grow quickly compared to the shear rate 
provided the vertical wavelength is much shorter than the transition wavelength Atran- Shear is more important near 
or longer than the transition wavelength; i.e. in the Parker regime, modes can be more affected by shear, unless they 
are associated with very short radial wavelengths. 

6. APPLICATIONS TO ASTROPHYSICAL ACCRETION DISKS 

In this section we investigate the relevance of photon bubble and Parker instabilities to the surface layers of high 
luminosity accretion disks, using recent stratified shearing box, r adiation-MHD simu latio ns as a guid e to ex pected local 
conditions. I n particular, we use data from simulations 0528a of Krolik et~al"1 (|2007| ) and lBlaes et all (|2007fl . and 1112a 
and 0519b of Hiro se et al.l ([2009). The simulations represent local patches of accretion disks at 30 (1112a and 0519b) 
and 150 (0528a) gravitational radii around a 6.62 solar mass black hole, and have time-averaged, volume-integrated 
radiation to gas pressure ratios of approximately 1 (0528a), 7 (1112a), and 70 (0519b). We time and horizontally 
average the fluid variables in each simulation in order to have an approximate vertically stratified equilibrium in which 
to explore the instabilities that we have calculated above. All time averages of the simulation data neglect the first 
ten orbits during which the magnetorotational instability was still developing. 

Figure [4] depicts the time and horizontally-averaged vertical profiles of radiation, gas and magnetic pressures in each 
of the three simulations. (The height z on the horizontal axis in these and other figures in this paper is in units of the 
scale height H of the initial condition of the simulations: approximately 3.1 x 10 6 cm for 0528a, 1.46 x 10 6 cm for 1112a 
and 4.37 x 10 6 cm for 0519b. This initial scale height is comparable to the actual total pressure scale height of the 
resulting time-averaged structure.) In the outer layers of the accretion disk in all three simulations, magnetic pressure 
greatly dominates gas pressure, and also generally dominates radiation pressure. As a result, the Parker instability 
generally controls the dynamics of these layers on the scales resolved by the simulations (jBlaes et al.ll2007l ). However, 
Figure 0] shows that at the highest levels of radiation to gas pressure ratio, radiation pressure is becoming comparable 
to magnetic pressure in supporting the disk outer layers. We therefore expect the photon bubble instability to be 
more prominent in higher radiation pressure systems as the transition wavenumber fctran decreases with decreasing 
magnetic to radiation pressure ratio, assuming comparable radiation and magnetic pressure scale heights. Photon 
bubbles, rather than Parker instability, may control the surface layer dynamics at high levels of radiation pressure 
support. 

Beyond pressure considerations, the medium also needs to be optically thick (but not so thick as to impede radiative 
diffusion) at a particular wavelength in order to drive the photon bubble instability at that wavelength. This is 
because diffusive radiat ive transport is essentia l for producing radiative amplification when the opacity is dominated 
by Thomson scattering (IBlaes fc Socr ates 2003). However, Figure [5] shows that the photon bubble turnover wavelength 
becomes optically thin to electron scattering beyond \z/H\ of approximately 2.1 and 2.7 for simulations 0519b and 
1112a, respectively. This means that the instability cannot exist at wavenumbers higher than fcx and thus will 
not achieve its maximum growth rate near the disk photosphere. These simple estimates mean that only longer 
wavelength photon bubble modes can grow in the upper layers of such accretion disks. Nevertheless, these slower, 
longer wavelength modes mi ght still develop in to shock trains: approximate nonlinear solutions with finite optical 
depths have been explored by iBegelmanl (|2006bO . The Parker instability, on the other hand, is not driven by radiative 



diffusion and hence is not subject to such optical depth restrictions. 

Incidentally, Figure [5] also shows that fcx > kp > k tran in all three simulations. With the exception of 0528a 
(see Figure [9]) where the gas pressure in the disk upper layer is not negligibly small, the growth rate dependence on 
wavenumber therefore resembles that depicted in Figure As shown in that figure, the Parker instability transitions 
into the photon bubble instability at short wavelengths before magnetic tension would introduce a cutoff. Note, 
however, that the photon bubble instability would never reach the asymptotic growth rate depicted in that figure 
because the turnover wavelength 2n/kT is optically thin. 

Figure |6] shows the same three characteristic wavelengths as a function of height in the simulations, now scaled by 
the azimuthal grid cell size. Also shown are the characteristic wav elengths that dem arcate the rapid (A < Ar) and 
slow (A > As) radiative diffusion regimes for acoustic perturbations (jBlaes et al.ll2007D : 



2-kc ( AE \ ( e 



3c{Kp \e -\- AE J \ct 

and 



(42) 



As = -A R . (43) 
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Fig. 4. — Horizontal and time averaged radiation (red), magnetic (blue) and gas (green) pressures as a function of height in local shearing 
box accretion disk simulations 0528a (top left), 1112a (top right) and 0519b (bottom). Vertical dashed lines indicate the locations of the 
Rosseland mean (mostly electron scattering) photospheres in the time-averaged structures. Note that the radiation to magnetic pressure 
ratio in the disk surface layers increases with increasing overall (volume integrated and time averaged) radiation to gas pressure ratio. 



Most of the results of this paper only apply in the rapid diffusion regime (A < Ar, i.e. below the green curves in 
Figure [5J We therefore see that Parker instabilities with small k x cannot exist in the rapid diffusion regime in the 
highest radiation to gas pressure ratio simulation 0519b, as they have already transitioned into photon bubbles. This 
reflects the fact that we noted above that magnetic pressure is not dominant over radiation pressure in the surface 
layers of this simulation. However, large k x modes (i.e. those with large radial shear between the field lines) will still 
be Parker-like in character even in the rapid radiative diffu sion regime, and the magnetic field and density structure 
in this simulation still shows evidence of Parker instability (jBlaes et al.H201lT ). 

Photon bubbles can reach their maximum growth rate in a narrow range of depths in simulations 0519b and 1112a, 
where radiative damping becomes small compared to the photon bubble asymptotic growth rate and th e turnover 
wave length is optically thick to scattering. Here the photon bubbles may evolve into nonlinear shock trains (|Begelmanl 
f200lb that can propagate outward into the photosphere and potentially produce observable signatures. However, 
Figure O shows that the turnover wavelength in 1112a and 0519b is at most a few times larger than the grid zone 
size in the region where &t is in the rapid diffusion limit, which means that the simulation cannot resolve the fastest 
growing photon bubble modes, even if they can physically exist in this region. 

Figure [7] shows the asymptotic, short wavelength, photon bubble growth rate (black curve), as well as the rapid 
diffusion Parker growth rate (blue curve), in the radiation dominated simulations 1112a and 0519b. The plots are 
restricted to a range of heights where radiative damping of photon bubbles is small enough so that they are unstable, 
and where the turnover wavelength is optically thick. Within this range, we only show the Parker growth rates where 
the cutoff wavelength Ap is in the rapid diffusion regime, where the analysis of this paper is valid. Photon bubbles 
grow faster than Parker in this range of heights, and so in principle could be important, but as we noted above, the 
simulations cannot resolve these fastest growing wavelengths. Note that the midplane regions, where most of the 
magnetorotational turbulence is acting, are not photon bubble unstable. 

Figure [8] shows the same information for simulation 0528a, and may explain why iBlaes et al.l (|2007f) did not observe 
the photon bubble instability in this P ra d ~ P gas simulation. There the transition wavelength in the upper layers 
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Fig. 5. — Optical depth of the characteristic wavelengths to electron scattering as a function of height for the horizontally and time- 
avcragcd structures in simulations 0528a (top left), 1112a (top right), and 0519b (bottom). Red, black and blue curves represent the photon 
bubble turnover wavelength, Parker cutoff wavelength and photon bubble/Parker transition wavelength, respectively. 

is moderately optically thick to electron scattering, which suggests that photon bubbles can physically exist in the 
simulated disk. However, Figure |S] shows that Parker is the dominant instability, at least for \z/H\ > 2.5, where Ap 
is in the rapid diffusion regime. Figure [9] shows the growth rate as a function of wavenumber at a representative 
height (z/H — 3.5) in the disk upper layer and further illustrates this conclusion. Closer to the midplane, radiative 
damping dominates over photon bubbles in the short wavelength limit. N ote that our plotted results are comparable 
to and consistent with the growth rates estimated by iBlaes et al.l (|2007[ ) at two specific epochs in this simulation. 
Those authors did not know how Parker and photon bubbles couple together, but we have demonstrated here that the 
maximum growth rates of the individual instabilities are in fact unaltered. 

7. CONCLUSIONS 

We expanded previous studies of photon bubble instabilities to include background magnetic field gradients, thereby 
introducing a coupling between photon bubble and Parker instabilities. Our first main result was obtained assuming 
negligible gas pressure in the equilibrium, which in principle allows the photon bubble and Parker instabilities to exist 
out to very short wavelengths. There we identified a finite transition wavenumber fct ra n between the photon bubble 
and Parker instabilities. The instability is Parker-like for k < fctran and photon bubble-like for fc > fctran) although 
strong horizontal shear between neighboring magnetic field lines (i.e. large k x ) pushes the Parker- like range to higher 
wavenumber. 

We then proceeded to a naive WKB study including finite gas pressure. We numerically found approximately the 
same transition wavenumber. Moreover, finite gas pressure introduces finite photon bubble turnover and Parker cutoff 
wavenumbers (kr and fcp, respectively). For fcp < fctran, we see a sharp growth rate decrease above k ~ kp, while the 
photon bubble growth rate becomes constant in k for k > kr- Once again, the growth rates of both instabilities agree 
with previous results in all wavenumber regimes. We also derived the scaling of the important wavenumbers fc-r, kp 
and fctran as functions of background parameters in equations (|38 j) -(|40 |) . Note that the results of this paper can only be 
applied to backgrounds where the characteristic wavenumbers fctran, fcp and fcx satisfy the WKB condition kH 3> 1. 

We also developed a WKB analysis of the coupled Parker/photon bubble problem including the effects of differential 
rotation. Photon bubbles are not significantly affected by shear and rotation provided their wavelengths are short 
enough. On the other hand, Parker modes can be affected by shear unless they have short radial wavelengths. This 
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Fig. 6. — The critical wavelengths Atran (blue), Ap (black) and At (red), normalized by the azimuthal grid cell size, as a function of height 



for the horizontally and time-averaged structures of simulations 0528a (top left, P ri 



3 ), 1112a (top right, P rad S 7P gas ) and 0519b 



(bottom, P ra d = 70P gas ). The green curves indicate the wavelength Ar below which radiative diffusion is rapid for acoustic perturbations, 
and the yellow curves indicate the wavelength Ag above which radiative diffusion is slow for acoustic perturbations. The gaps in the black 
curves indicate regions where the horizontal and time averaged magnetic pressure increases vertically outward. Our WKB analysis is only 
valid below the horizontal dashed lines, which indicate the value of 2nH in each of the simulations. 





Fig. 7. — Maximum rapid diffusion Parker gr owth rate (blue) [f r om eq uation dA13M and the asymptotic short-wavelength photon bubble 
growth rate [black, based on equation (93) of Blacs & Socrates (2003)] as a function of height from horizontal and time averaged data 
above the midplane in simulations 1112a (left) and 0519b (right). The plotted growth rates are scaled with the local orbital period for each 
simulation. The wave vector orientation is k = (0, cos(-7r/4), sin(-7r/4)). At heights below those shown in the plots, photon bubble growth 
is suppressed by radiative damping, while at heights above those shown, the photon bubble turnover wavelength is optically thin, so that 
the short wavelength asymptotic growth rates cannot be achieved. Note that the plotted Parker growth rate is only valid within regions 
where the cutoff wavelength kp is in the rapid diffusion regime, which in the case of 0519b is above the plotted range. 




Fig. 9. — Instability growth rate as a function of wavenumber for parameters taken from the horizontal and time averaged simulation 
0528a results at z/H = 3.5. We see that the maximum Parker instability growth rate is indeed higher than the photon bubble asymptotic 
growth rate in this case, in agreement with Figure [8] 
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is equivalent to the condition of strong horizontal shear between neighboring held lines that expands the Parker-like 
range to higher wavenumber. 

We applied our results to analyze radiation MHD shearing box accretion disk simulations 0528a, 1112a and 0519b 
with volume integrated radiation to gas pressure ratios of 1, 7 and 70, respectively. We found that the asymptotically 
growing short wavelength photon bubble instability is not likely to exist in the upper layers of both 1112a and 0519b 
because the turnover wavelength becomes optically thin to electron scattering in the photosphere. For 0528a, the 
undulatory Parker instability in the rapid diffusion regime dominates over photon bubbles in the upper layers. On the 
other hand, radiative damping destroys the photon bubbles near the midplane in all three simulations, and photon 
bubbles may exist at the asymptotic growth rate only in a narrow range of heights. Even there it is not clear whether 
the turbulent conditions would allow the photon bubbles to grow unimpeded. Last but not least, the grid cell sizes of 
current radiation pressure dominated simulations are inadequate to fully resolve the turnover wavelength of the photon 
bubble instability, so that the fastest growing wavelengths cannot manifest themselves. Higher resolution radiation 
pressure dominated simulations may reveal this instability. 

While accretion disk applications have been our focus, much of the analysis of this paper applies to static media, 
and magnetize d stellar enve l opes are clearly another area where this physics would be relevant. Generalizing the 
simulations by iTurner et al.l (j2005| ) of magnetized, vertically stratified static equilibria to include magnetic pressure 
gradients would be interesting in order to explore how the photon bubble and Parker instabilities interact in the 
nonlinear regi me. Even if short wa velength photon bubbles have faster linear growth rates than the longer wavelength 
Parker modes. [Turner et al.l (|2005l ) have shown that the resulting nonlinear shock trains merge, causing the distances 
between adjacent shocks and the density contrast to grow until the magnetic field buckles (or until only one wavelength 
fills the simulation domain) . This field line buckling may resemble the nonlinear development of the undulatory Parker 
instability, and it may be that both instabilities ultimately lead to a similar nonlinear outcome, at least in some 
regimes. 

Finally, we stress an important caveat to our work: the photon bubble-Parker transition regime of our numerical 
finite gas pressure results do not readily apply when c\ > i>a, c r or c r > i>a, due to limitations of the WKB method 
discussed in section 4. Such regimes can exist in the portion of the accretion dis k within and just o utside the MRI 
unstable region, where the gravitational potential energy is dissipated into heat ()Hirose et al.N200"9l ). The behavior 
of the Parker-photon bubble transition in this region may be important for understanding how accretion power is 
transported outward. This is a difficult problem that may deserve further future work. 

We owe a debt of gratitude to Ellen Zweibel for illuminating many aspects of the physics of the Parker instability. 
We also thank the referee for constructive criticism and for suggesting that we look more carefully at the issues of 
rotation and shear when considering accretion disk applications. We have benefited from helpful discussions with L. 
Bildsten, M. Block, K. Choiu, E. Gallo, J. Jacob, E. Rykoff, N. Turner, E. Newton and D. Harsono. This work was 
supported in part by NSF grants AST-0307657 and AST-0707624. 



APPENDIX 



MAGNETIC BUOYANCY INSTABILITIES IN THE ABSENCE OF PHOTON BUBBLES 



In order to better understand the coupled photon bubble-Parker problem, it is useful to derive the properties 
of pure Parker modes in the absence of photon bubble instabilities. We are most interested in a medium where 
radiation pressure dominates gas pressure, and in short wavelengths where radiative diffusion is rapid. The latter 
condition implies that temperature perturbations are very small, and if we take them literally to be zero, then the 
flux perturbations that drive photon bubble instabilities in this regime will vanish identically, leaving only pure Parker 
modes. 

Parker modes with rapid radiative diffusion were first studied by Gilmanj (|1970f ). who argued that it was a good 
first order approximation to replace the perturbed energy equation (] 14|) with the statement that the temperature 
perturbation ST vanishes. 

The linearized perturbation equations that we presented in subsection 2.2 can be combined with the perturbed flux- 
freezing equation to give the total (gas plus radiation plus magnetic) pressure perturbation in two radiative diffusion 
limits. The first is the adiabatic limit of infinitely slow diffusion (k —¥ oo, F —¥ 0, kF finite), giving 



SPu 





k 2 y v\ 








UJ 2 




\c 2 N 2 


UJ 2 J 


9 




c(N 2 



o d 
v 2 A —lnB 
A dz 



2 d (B 

dz \ p 



(Al) 



where £ z = j 5v z dt = iSv z /uj is the vertical component of the Lagrangian displacement vector. The second is the 
limit of sufficiently short wavelengths that diffusion is rapid enough (k — > 0) to guarantee that the perturbations are 
isothermal (8T — > 0). This gives 
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Equations (|A1|) and (|A2[) contain the essential physics of the Parker instability in both of these limits. Buoyancy 
is maximized for a vertically displaced fluid element if it is in pressure equilibrium with its surroundings, i.e. SPtot — 
0. Solving the resulting equation for the Eulerian density perturbation associated with an upward fluid element 
displacement (£ z > 0) then gives instability if 8p < 0. There are two possibilities depending on k y . For k y — > 0, the 
first £z term on the right hand sides of equations (IA1|) and (|A2|) dominates. In this case, the perturbation involves 
a straight bundle of field lines, and the resulting mode is called an interchange mode. On the other hand, for k y 
sufficiently large, the second £ 2 term dominates, and the perturbation involves a bending of the field lines, resulting 
in the undulatory mode. 

In t he adiabatic limit, the instability criteria for these two modes (jNewco mb 196H lAchesorill979tlChristodoulou et al.l 
|2003f ) are therefore 
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In the rapid diffusion limit, on the other hand, the instability criteria may be written as (Gi lmanl 119701 : lAchesonl 
1979) 
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These are easier to satisfy than the criteria (IA3[) because pressure and temperature equilibrium between the pertur- 
bations and their surroundings makes the fluid hydrodynamically neutrally buoyant. 

We consider only the rapid diffusion limit from now on. Using equation (|A2[) . the linearized continuity and horizontal 
momentum equations can be combined to express the density and total pressure perturbations entirely in terms of the 
vertical velocity perturbation: 
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Here k ± = (k 2 x + k 2 ) 1 / 2 is the magnitude of the horizontal wavenumber. 
These expressions can then be combined with the linearized vertical momentum equation, 
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to give an exact linear ordinary differential equati on for the verti cal velocity perturbation. In his analysis of the Parker 
instability in the limit of rapid radiative diffusion, iGilmanl (|1970f) derived a simplified form of this differential equation 



by first assuming k x 
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oo with uj and k y remaining finite. From the x-component of the linearized momentum 
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this guarantees that SPtot 0, giving a maximum buoyant response and therefore a maximum Parker growth rate. 
However, here we will allow for the possibility that k y /k x is not necessarily small in magnitude as k x gets large. This 
is because the Parker cutoff wavenumber fcp = [g/(2c?_ff mag )] 1 / 2 can be very large in a radiation dominated plasma, 
in which the isothermal sound speed in the gas alone c; can be very small. (In contrast, in the adiabatic limit, the 
Parker cutoff wavenumber is comparable to the inverse scale height of the background medium.) Because a plasma 
whose thermal pressure is dominated by radiation is necessarily supported against gravity by radiation pressure and/or 
magnetic pressure gradients, the Parker cutoff wavenumber in the rapid diffusion limit can be much larger than the 
inverse scale height of the background. It therefore makes sense to consider the possibility of very large, as well as 
very small, k y . 
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We therefore take \uj 2 \ <C kj_v A , but not necessarily much smaller than k 2 v\. The first derivative term in equation 
(|A5|) is higher order than the other terms for short wavelengths, so the density perturbation becomes 
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For the total pressure perturbation in equation (|A6j) . the first derivative term dominates, and 
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Substituting into the vertical momentum equation, employing the short wavelength WKB approximation d/dz 
ik z , and simplifying, we finally obtain the dispersion relation 
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where k = (fcjj_ + k 2 ^) 1 ! 2 is the total wavenumber magnitude of the perturbation. Apart from the generali z ation that 
k y and k z no longer need be consider ed sm all compared to k x , this is identical to equation (14) of Gi lmanl (|1970f) . 

If either of the instability criteria (IA4|) arc satisfied, then all nonzero wavenumbers k y less than k±kp/k will be 
unstable, where, again, the characteristic Parker cut off wavenumber is defined by 



kp 



2r 2 H 



The maximum undulatory instability growth rate is given by 
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Assuming that the gas sound speed is much less than the Alfven speed (the regime of interest for this paper), 
equations (|A13|) and (|A14I) become 
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THE ZERO GAS PRESSURE LIMIT WITH ROTATION AND SHEAR 

Here we generalize the zero gas pressure limit analysis of section [3] to include the effects of rotation and shear on the 
coupled Parker/photon bubble instability problem in a differentially rotating accretion disk. Because we are interested 
in short wavelength perturbations, we focus on a local comoving patch of the accretion disk, and employ the shearing 
box approximation (|Goldreich fc Lvnden-Belll 119651 : iHawlev. Gammie fc Balbusl I1995D . The fluid equations ((U)-© 
remain the same, except for two modifications. 

The first is the addition of the terms — 2pttz x v + 2qpQ 2 xyL to the right hand side of the momentum equation @, 
which represent the combined effects of Coriolis, centrifugal, and radial gravitational forces on the flow. Here fi is the 
angular velocity of the local patch of disk and q is a shear parameter, representing a variation of angular velocity with 
radius r of ft oc r~ q (q — 3/2 for Keplerian shear). We employ a Cartesian basis (x, y, z), with x in the local radial 
direction, y in the local azimuthal direction, and z in the vertical direction just as in the static equilibrium studied in 
the main body of the paper. We continue to write the vertical gravitational acceleration as g = — gz, with g = tt 2 z, 
the equilibrium radiation flux as F = Fz, and the equilibrium magnetic field as B = B(z)y, i.e. a purely azimuthal 
magnetic field. The equilibrium state continues to be described by equations but in addition there is now an 

azimuthal flow velocity of v = — qVLxy. 

The second modification concerns horizontal boundary conditions. In the static analysis, we considered a horizontally 
homogeneous equilibrium, and perturbations that were periodic in the horizontal direction (cx expf^fc^a; + k y y)]). 
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The shearing box has all fluid va riables subject to shearing periodic boundary conditions in the radial direction 
(|Hawlev. Gammie~fc Balbus"l ll995f) . The standard way to treat this is to make a coordinate transformation: x' = x, 
y' = y + qflxt, z 1 = z, and t' = t. Our perturbations can then be taken to have a horizontal dependence of the form 
exp[i(k x (t')x + k y y)\, where the wavenumber in the x-direction is time-dependent: 

k x (t') = k x o + qnt'ky (Bl) 

Here k x o and k y are independent of time. Horizontal derivatives acting on perturbations can therefore be written as 
d/dx —> ik x (t') and d/dy — > ik y . We also adopt the Lagrangian time derivative operator 

d d d 

J?"S-«*5- (B2) 

The explicit time-dependence of the radial wavenumber implies that a Lagrangian treatment of the perturbation 
equations is much easier than an Eulerian treatment. We therefore work in terms of the Lagrangian displacement 
vector which to linear order is related to the Eulerian velocity perturbation by 

<5v = -± + gfi^y. (B3) 

This allows us to immediately solve the continuity and flux-freezing equations to exp ress the Eulerian density and 
magnetic field perturbations in terms of the Lagrangian displacement (lNewcomblll962f ). To linear order, 

S P = -ipk x £ x ~ ipkyiy - p^p - £ z ^p (B4) 
5B x =ik y B£ x , (B5) 

5B y = -Z z ^--B^-iBk x Z X) (B6) 
dz dz 

and 

6B Z = ik y B£ z . (B7) 
Eliminating the density and magnetic field perturbations, the linearized momentum equations are then 

= -ijSP - (kl + fc>^ + 2 ^ + 2 "§7 - + ^% (B8) 

ih 5P 2i & (B9) 



dt' 2 p dt' 2H> 



mag 



and 



d% _ 1 dSP - a vj 3vj d£ z 2 d% .k x v% 

dt' 2 ~ p dz > VA ^ Z + {H' m ^ z 2H mas dz +VA dz 2 l H m J x 



+ ik x v A~Q^ + i9 k x£x + igky^y + g - (BIO) 

The thermal pressure perturbation can be derived from the linearized energy equation. Assuming negligible gas 
pressure, and eliminating the flux perturbation using the radiation diffusion equation, this is 

d5P dE A r dSp c dSP c ,,o , 2 , r c d 2 SP 

Thus far, all of our equations are exact in the zero gas pressure limit. We now employ the WKB approximation in 
z, assuming that all perturbations have a ^-dependence of exp(ifc z z) and take all wavenumber components to be large. 
Because we are looking for instabilities with growth rates that increase with wavenumber no faster than fc 1 / 2 , the time 
derivative term in the energy equation is small. Physically, radiative diffusion is extremely fast at short wavelengths 
compared to our instability growth times, i.e. we are in the rapid diffusion limit. Eliminating the time derivative 
here eliminates the damped diffusion mode of equation (|19[) . The remaining terms with the largest dependence on 
wavenumber then allow us to approximate the thermal pressure perturbation as 



6P^^ 
ck 2 



Fk z {k x ^ x + ky^y + k z t; z ) — (ik x + ik y ^ + ik z + ik y qfl^ x 



(B12) 



The momentum equations can also be simplified by first noting that the second order time derivative terms on 
the left hand sides of equations (IB8|) and (IB10[) are negligible given the ~ k 2 magnetic force terms on the right 
hand sides. Physically, eliminating these time derivatives eliminates the fast and Alfven modes of equations (|20[) and 
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(p?T]) . Equations (|B8I) and (IB10[) can then be used to eliminate £ x and £ z from equation (|B9[) . The algebra is greatly 
simplified if one first recognizes that in this limit of negligible gas pressure, slow modes have fluid displacements that 
are predominately along the magnetic field (in the y-direction) . In fact, £ x and £ z are one power of k smaller than 
Keeping only the lowest order terms as k — ¥ oo, and assuming that magnetic gradients dominate radiation pressure 
gradients in supporting the equilibrium (so that k tran H mag S> 1 and the Parker regime generally exists in the WKB 
limit for all wavenumber orientations), we finally obtain the following differential equation: 

Si' 2 + fc 2 3c dt> AHl ag fc 2 A: 2 c ~ 1 j 

If there were no shear (q = 0), so that the radial wavenumber fc^ was independent of time, we could could assume a 
complex exponential time-dependence £ y cx exp(— iut') and recover the static equilibrium dispersion relation (|27[). 

The Coriolis accelerations are simply too small to have affected the modes at short wavelen gths . To orde r of 
magnitude, assuming a typical Parker growth rate ~ VA/H mag , the Coriolis terms in equations (|B8[) and (IB9|) are 
negligible provided 57 <C kv&- But for a magnetically supported equilibrium, ~ v&_/H mag , so this is just equivalent 
to the WKB condition kH mag >• 1. In more detail, we require \k z \H mag ^> k \k v \/(\k x \k ±), so that large radial 
wavenumbers better ensure immunity of the Parker instability from Coriolis effects (|Shulll974l) . Provided k tran H mag 3> 
1, the photon bubble regime will be even less affected by Coriolis forces. 
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